A deep XMM-Newton survey of M 31* 
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Abstract. The deep homogeneous survey of the large Local- 
Group spiral galaxy M 31 is a milestone project for X-ray as- 
tronomy, as it allows a detailed X-ray inventory of an archety- 
pal low-star-formation-rate galaxy like our own. We present 
first results of the survey, which covers the entire D25 ellipse. 
Information from different X-ray energy bands are combined 
in an X-ray colour image of M 31. In the first 15 observations 
we found about 1000 sources, the full survey will yield about 
2000 X-ray sources. 

Sources will be classified using hardness ratios, extent, high 
quality spectra and time variability. In addition the sources will 
be correlated with catalogues in optical, infra-red and radio 
wavelengths. 

Our goal is to study M 3 1 X-ray binaries and globular cluster 
sources, supersoft sources, supernova remnants and the hot in- 
terstellar medium and separate them from foreground stars and 
background objects. 



1. Introduction 

The goal of this survey is to improve our understanding of the 
different source classes in M 31 and possible connections to 
the general structure of the galaxy. Additionally the interstellar 
medium (ISM) will be studied. The results achieved by this sur- 
vey will help us to better understand the source classes of the 
Milky Way and of more distant galaxies . Due to its relatively 



short distance (780 kpc, lHollandlll998l: IStanek & Garnavich 
19981) and its moderate foreground absorption (A^h= 7x10^° 



cm 



Stark et al.lll992h . the bright, massive spiral (SA(s)b) 



galaxy M 3 1 is an ideal target for such a study. 

Early observations in the X-ray wavelength regime were 



Send offprint requests to: H. Stiele, email: hstiele@mpe.mpg.de 

* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 



perfo r med by the E i nstein observa tory (van Sp evbroeck et al . 
I979I: ICollura et al.1 Il990l: iTrinchie ri & Fabbiano" [l99lb and 
ROSAT jPrimini et al.l Il993t ISupper et al.l Il997[l200ll) . In 
the last few years many observations with Chandra, XMM- 
Newton, Swift and other X-ray satellites covered mainly the 
central region of M 31. Individual references will be given 
in Section |4] where the different source classes are discussed. 
Cha ndra resuhs are discussed in a different paper of this meet- 
ing ( IWilliamsll2007l) . 



2. Observations and data 

The survey consists of 15 individual EPIC observations. They 
are placed in such a way that the whole D25 ellipse is covered 
(see Fig. [1]). The observations were taken in XMM-Newton 
AO 5, between June 2006 and January 2007. 
For each observation a high energy (7 to 15 keV) background 
light curve was produced. We used these light curves to se- 
lect time intervals with low background, the so-called good 
time intervals (GTIs). As many of the observations were af- 
fected by strong background flares (mostly at the beginning 
or end of the observation), the net exposure useful for our 
analysis is strongly reduced. The sum of GTIs after high en- 
ergy background screening ranges from 6 to about 55 ks. The 
observations taken during the summer visibility window of 
M 3 1 were affected more strongly by background radiation than 
those taken during the winter window. The most affected ob- 
servations were reapproved in XMM-Newton AO 6. A severe 
background screening is especially important for the detection 
of faint point sources, investigation of diffuse emission and for 
colour images. We also want to mention, that for some of the 
observations the soft energy (< 2 keV) background varied and 
an additional screening was necessary. 

We used five energy bands: (0.2-0.5) keV, (0.5-1.0) keV, (1.0- 
2.0) keV, (2.0-4.5) keV, and (4.5-12) keV, to create images, 
background images and exposure maps for PN, MOS 1 and 
MOS 2 and masked them for acceptable detector area. For PN 
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Fig. 1. This image shows the pointings of the 15 observations 
of the large survey overlaid on a deep optical image of M 31. 
The region covered by a single observation is marked with a 
bold circle with a radius of 14 arcmin. 



Fig. 2. First combined EPIC PN, MOS 1 and MOS 2 RGB im- 
age of the new and archival data. For more details see Sec. [3] 



the background maps contain the contribution from the "out of 
time (OOT)" events. 

Source detection is done on all 15 images simultaneously, 
using the XMM-Newton Science Analysis System (SAS) 
tasks eboxdetect and e mldetect. Details can be found 
e.g. in IStiele et al.l (l2007al) . We detected about 1000 sources 



in the A05 o bservations. Together with the sources found 
by lPietsch et a l. (2005b) and the reapproved observations, the 
XMM-Newton source catalogue will contain about 2000 indi- 
vidual sources in the field of M 3 1 . 

3. Colour image 

Figure |2] shows a first combined exposure corrected EPIC PN, 
MOS 1 and MOS 2 RGB image of the new and archival data. 
The colours represent the X-ray energies as follows: red 0.2 - 1 
keV, green 1-2 keV and blue 2-12 keV. The optical extent 
of M 31 is indicated by the D25 ellipse and the boundary of 
the observed field is given by the green contour The image is 
smoothed with a Gaussian of 20" FWHM. Observation S3 was 
removed completely due to strong background. In some obser- 
vations individual noisy MOS 1 and MOS 2 CCDs are omitted. 
The colour of the sources reflects their class. Supersoft sources 
(SSSs) appear in red. Thermal supernova remnants (SNRs) and 
foreground stars (fg stars) are red to yellow. "Hard" sources 



(background objects (AGNs), X-ray binaries (XRBs) or Crab- 
Uke SNRs) are blue. 

4. Source classes 

We want to classify and if possible identify the sources of our 
survey. In this way we will separate sources belonging to M 3 1 
from foreground stars and background objects. In the following 
we will discuss the different methods, which help us to differ- 
entiate the source populations. 

4.1. Hardness ratios 

The most general method, which can be applied to all sources 
and which is based on their X-ray spectral properties, is the 
analysis of hardness ratios. The hardness ratios and errors are 
defined as 



HRi = 3^ and EHRi = 2 — ir-^; ,(1) 



Bi+i + Bj 



for / = 1 to 4, where B,- and EBj denotes count rates and cor- 
responding errors in band / as defined above. A hardness ratio 
can be recognised as an X-ray colour 

Figure[3]shows HR2 plotted versus HRI for sources of the cat- 
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Fig. 3. Hardness ratio plot HR2 ve rsus HRl for the sources of 
the catalogue o fiPietsch et al.l(l2005bi) . Marked are the three re- 
gions, which separate clearly from each other. Source classifi- 
cation is indicated: Foreground stars and candidates are marked 
as big and small stars, AGN candidates as small crosses, SSS 
candidates as triangles, SNR and candidates as big and small 
hexagons, GlCs and XRBs and candidates as big and small 
squares. 



alogue of iPietsch et al.l (l2005bl) . The sources group in three re- 
gions, which separate clearly. In the lower left corner SSSs can 
be found (HRl < and HR2 ~ -1). The major class of SSSs 
in M 3 1 are optical novae. This was found by iPietsch et al. 
(l2005al) correlating SSSs detected with ROSAT, Chandra and 
XMM-Newton with optical nova catalogues from the litera- 
ture and from their own observation s. The project is con tinued 
and recent results are published in iPietsch et al.l ( l2007b . Two 
SSSs which are not associated with optical novae are known 
to show periodicities of 865 s and 217 s, respectively. They 
are suggested to be spinning white dwarfs in a post-nov a state 



(lOsborne et al ] l200lt iTrudolvubov & Priedhorskyl2007h . 
SNRs and fg stars are located in the lower right corner of that 
hardness ratio plot (HRl > -0.1 and HR2 < -0.2). All "hard" 
sources (AGNs, XRBs and Crab-like SNRs) can be found in 
the upper right. 

4.2. Time variability 

To further distinguish these source classes, an investigation of 
their X-ray time variability is important. SNRs do not show 
time variability of their overall intensity on time scales of years. 
In this respect SNRs and fg stars can be separated which is not 
possible from HRs alone. A fg star classification can be sup- 
ported by an optical counterpart or if the star also shows vari- 
ability on short timescales. A reason for this variability can be 
flares (see Fig.|4]and below). 

With the help of time variability one can also separate XRBs 



from AGNs and Crab-like SNRs. Only XRBs can show strong 
time variability. Among the AGNs we know only a handful in 
the whole sky which show variability comparable to that of 
XRBs. So it is very unlikely that a strongly (Fyai > 10) vari- 
able source in the field of M 31 would be an AGN. A detailed 



study o n the central region of M 31 is presented by lStiele et al 
(l2007ah . 



In addition XRBs can be transient sources. T here are many 



studies o f transient sources in M 31 Uke e. g. IWilUams et al 



('2006b'), 'Trudolvu bovetaD (l2006h . IWilhams et al.l (l2005h . 



iWilliams et al.i (i2006al) . A transient which we examined within 
our survey is discussed in Section |5] It should be mentioned 
that up to now among XRBs only low mass systems are known 
inM31. 

4.3. Cross correlations 

So far only X-ray properties of the sources have been used for 
source classification. Additional information can be achieved 
using cross correlations with catalogues in the radio, infra- 
red and optical wavelengths. X-ray sources which correlate 
with optical globular clusters are most likely X-ray bina- 
ries. Bright XRBs in globular clusters were investigated by 
Di Stefano et al.l (|2002|) . Optical images can be helpful to dis- 
tinguish between fg stars and SNRs. Correlations with radio 
counterparts can support AGN and SNR classifications. For the 
classification of AGNs also optical spectra can be helpful. In 
Chandra obser vations SNRs can b e found as spatially resolved 
X-ray sources dKong et al.ll2003b . So the extent of an X-ray 
source can be an additional criterion for a SNR classification. 



4.4. Additional methods for bright sources 

We examined 1000 s light curves of all bright objects in our 
survey for flares. An example is given in Figure|4] If in addition 
an optical counterpart with low ratio of the X-ray to the optical 
(V band) flux (log(/x//opt) < -1.0) is found, the source must 
be a fg star To clearly see the rise and decline of the flare it is 
necessary to have a long and continuous observation. 
The classification of brighter sources can be improved using 
X-ray spectra and short term time variability. 



5. A new X-ray transient 

The transient (CXOM31 J004 059.2-H41 1551) was detected 
with Chandra on 2007 July 5 (iCalache et al .1120071) . The po- 
sition of the transient nearly coincides with the centre of the 
SN 1 pointing of our survey. During a 22 ks XMM-Newton tar- 
get of opportunity observation (TOO) on 2007 July 25, about 
20 days after the Chandra detection, the source was still bright. 
We fitted the 0.2 - 7.0 keV EPIC PN spectrum (see FiglD by an 
absorbed disk blackbody model with best fit parameters A^h = 
2.1+0.2xl02'atoms/cm2 and Tin = 0.51 +0.02 keV. The unab- 
sorbed luminosity (0.5 - 10 keV) is ~ 9.1 x lO^^erg/s. The 
spectral parameters and luminosity are in agreement with the 
Chandra values. A FFT period search did not reveal any sig- 
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Fig. 4. Combined EPIC PN and MOS 1000 s light curve of a 
bright source showing a flare. Corr esponding counte rparts are 
found in the Local Group Survey (iMassev et al .l2006h . 2MASS 
and USNO-Bl catalogues. Thus it can be cla ssified as fg star. 
This source was also detec ted with Chandra iPi Stefano et al 
2004 nl-15) and ROSAT dSupper et al.|[l99^ 
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Fig. 5. 0.2 - 7.0 keV EPIC PN spectrum of the transient 
CXOM31 J004059.2+411551, fitted by an absorbed disk 
blackbody model. Best fit values given in Sec. ID 



nificant periodicities in the 0.3 to 2000 s range. We did not find 
any evidence of optical counterparts for the transient, check- 
ing the XMM-Newton optical monitor UVWl, UVM2 images, 
taken during the T OO observation, ari d the Local Group Sur- 
vey M 31 images dMassev et al. I l2006h . The X-ray parameters 
and the lack of an optical counterpart are consistent with th is 
source being a black hole X-ray transient (iStiele et al ] l2007bh . 



6. Summary 

The deep XMM-Newton survey of M 31 consists of fifteen 
pointing directions, which are placed in such a way that to- 
gether with the observations in the XMM-Newton archive the 
whole D25 ellipse of the galaxy is covered. The observations 
were taken in AOS between June 2006 and January 2007. Some 



of the observations affected by high background were reap- 
proved in A06. The full survey will yield about 2000 sources in 
the field of M 3 1 . This survey will separate sources belonging 
to M 3 1 (supersoft sources, quasisoft sources, supernova rem- 
nants and X-ray binaries) from background objects, especially 
AGNs, and foreground stars. 

We produced an RGB colour image in which SSSs appear as 
red, fg stars and SNRs as red to yellow sources. All "hard" 
sources are blue. A more qualitative form of the selection cri- 
terion, which forms the basis of this colour selection, are hard- 
ness ratios. In addition X-ray time variability can be helpful to 
distinguish between different source classes. Further informa- 
tion is achieved by cross correlations with catalogues in other 
wavelengths. For bright sources X-ray light curves. X-ray spec- 
tra and short term time variability can be examined. 
The next steps will be an improvement of the source positions 
and the preparation of a source catalogue. 
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